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We study the spin-dependent electronic transport through a one-dimensional ballistic quantum 
wire in the presence of Rashba spin-orbit interaction. In particular, we consider the effect of the 
spin-orbit interaction resulting from the lateral confinement of the two-dimensional electron gas 
to the one- dimensional wire geometry. We generalize a situation suggested earlier [P. Streda and 
P. Seba, Phys. Rev. Lett. 90, 256601 (2003)] which allows for spin-polarized electron transport. 
As a result of the lateral confinement, the spin is rotated out of the plane of the two-dimensional 
system. We furthermore investigate the spin-dependent transmission and the polarization of an 
electron current at a potential barrier. Finally, we construct a lattice model which shows similar 
low-energy physics. In the future, this lattice model will allow us to study how the electron-electron 
interaction affects the transport properties of the present setup. 

PACS numbers: 72.25.Dc, 71.70.Ej, 72.25.Mk 
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I. INTRODUCTION 

Spin-orbit coupling is a relativistic effect of order 
0(v 2 /c 2 ), where v is the electron velocity, which follows 
directly from the Dirac equation. It is described by the 
Hamiltonian (for V x E = 0) 



Hso = 



eh 



Am 2 c 2 



E x 



(1) 



where the electric field E = —W/e (e < is the elec- 
tron charge) is the gradient of the ambient potential. 
In the following, the correction — eA/c to the canonical 
momentum is abandoned. In order to confine electrons 
to nanostructure devices, sharp potentials are necessary, 
which lead to nonnegligible spin-orbit interaction (SOI), 
especially in systems with structural inversion asymme- 
try like e.g. semiconductor heterostructures. This effect 
can be used to achieve control over the electron spin 
and leads to spin-dependent transport properties, such 
as spin-polarized currents, even in systems without fer- 
romagnetic leads. 

The emerging field of spintronics might result in an ex- 
tensive use of the spin degree of freedom for information 
processings^ In a two-dimensional electron gas (2DEG) 
obtained by a strong confinement in z-direction, the SOI 
is usually described by the so-called Rashba term 



Hr = —a z (cr x py - o y p x ) 
m 



(2) 



contributing to the Hamiltonian of the electron system.—^ 
Here the components of the electron momentum oper- 
ator are denoted by Pi, the Pauli matrices by <Tj, and 
a z oc E z is the SOI coupling coefficient^ set by the con- 
fining electric field. As discussed by Datta and Das^ a 
further confinement of the 2DEG to a wire geometry al- 
lows for a particular control over the spin, if a z or the 
length of the wire are varied. This insight led to exten- 
sive studies on the transport properties of noninteracting 
electrons in quasi one-dimensional (ID) quantum wires 



with SOI&l&SiiOiiLi^ In p art i cu i arj tne effect of sub- 
band mixingiiiSiii and a magnetic field perpendicular to 
the plane of the underlying 2DEG 12 was investigated. 

A very promising candidate for a system to experimen- 
tally produce spin-polarized currents using SOI is the 
setup suggested by Streda and Seba where the magnetic 
field points in the wire direction and an additional po- 
tential step is placed in the quantum wire.— ^ It is as- 
sumed that due to the large energy level spacing only 
the lowest subband of the quantum wire is occupied and 
subband mixing can be neglected. Restricting the con- 
siderations to this subband, one does not have to in- 
clude explicitly the potential confining the electrons to 
the wire. Furthermore, the strong lateral confinement 
allows to take into account only the momentum in the 
wire direction, p x = p, p y = p z = in Eq. @. The en- 
ergy dispersion of the ID electron gas £o(k) = h 2 k 2 / (2m), 
where k — k x , is split by the Rashba term Eq. ^ into 
two branches e^ s \k) 
s = ± and E ctz = 1 
fourfold degenerate with two left and two right moving 
states. The spin expectation values are (a y )k, s = s and 
{ox)k,s = {oz)k,s = 0, independent of In presence of an 
external magnetic field (parallel to the wire), described 
by a Zeeman term 



: h 2 (k + sa z ) 2 /(2m) - E az , with 
a 2 /(2m). The eigenenergies are 



H z = ezc.r/2 , 



(3) 



an "energy gap" of size ez opens up at k = [see Fig.Q]a)] 
and states within this "gap" are only twofold degenerate 
(one left and one right moving state). A potential step 
can then be used to generate a tunable spin polarization, 
in mainly the y-direction, of the linear response current. 
In order to achieve this, the height of the step Vo > 
for wire positions x < has to be chosen such that the 
energy falls into the "gap" region, while for the potential 
free part a; > it lies sufficiently above the "gap" [see 
Fig. [1] a)] . As an additional effect of the magnetic field, 
the spin expectation value is rotated gradually from the 
±y-direction into the ±x-direction when \k\ — > 0, while 
(cr z ) k s remains zero. Depending on the chosen param- 
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FIG. 1: (Color online) a) A potential step of height Vb and 
b) a potential barrier of height Vb and width 2x c . The cor- 
responding dispersions in the different regions are sketched 
(solid line: s = +, dashed line: s = — ). 



study the transmission current and the spin polarization 
at a potential barrier and discuss the interplay of a y and 
a z . In addition, we present a lattice model which in an 
appropriate parameter regime shows the same physics 
as the continuum model. This model will allow us to 
study the effect of the electron-electron interaction on 
the spin polarization in a forthcoming publication 1 ^ us- 
ing the functional renormalization group method. 14 



II. CONTINUUM MODEL 



eters, this leads to a small x-component of the ground 
state magnetization, whereas the y- and z-components 
are exactly zero as will be explained below. 

We here generalize the situation studied in Ref. H in 
several ways. We first study how the above scenario is 
modified in the presence of an additional Rashba term 

H' R = — a y a z p x (4) 
m 

resulting from the confinement of the 2DEG to the wire 
geometry, a term which so far was mainly ignored. As we 
also focus on the lowest subband and do not study sub- 
band mixing, the exact shape of the potential confining 
the electrons to the wire is not important. As its main 
effect, H' R will lead to nonvanishing spin expectation val- 
ues {(J z )k,s and thus a spin polarization component per- 
pendicular to the plain of the underlying 2DEG. We also 



The model we consider is given by the Hamiltonian 

pi ha z ha v eh _ 

H = -a vPx + &zPx ~ 7; a- B 5 

2m m m Zmc 

We slightly generalized the situation discussed above and 
allow for a Zeeman term with a magnetic field B = 
_B(sin 9 cos ip, sin 9 sin (p, cos 9) pointing in arbitrary direc- 
tion. The normalized eigenstates with quantum numbers 
k and s = ± are given by the product of a plane wave 
(in x-direction) and a two-component spinor 

Applying the Hamiltonian Eq. |5]) to this ansatz we ob- 
tain 



k 2 + 2a y k + 2fcfcos0- 
2ika z +2kle i * sin 9, 



2ika z + 2k%e~ i v sin6» 
fc 2 — 2a v k — 2k% cos 9 — e 




= 0, 



(7) 



with e = 2mE/h 2 , a y — eE y /(4mc 2 ), a z = eE z / (4mc 2 ), and k\ = —eB/(2Hc). Note that a y ,a z < in our notation 
due to the negative electron charge. One obtains the eigenenergy (divided by h 2 /2m) 



5^(k) = k 2 + 2ssgn(k - ko)y/C(k), 



(8) 



with C(k) — (a 2 + a 2 )k 2 + 2k 2 z k(a y cos9 — a z sin 9 sin <p) + k\ and fco = — k\ (a y cos 9 — a z sin 9 sin ip) /(a y + a 2 ) being 



the 



number at which the "energy gap" becomes smallest [see Fig. [2]. The corresponding eigenfunctions are 



» 



(9) 



with 



-ia z 



-up , 



Oiyk 



+ k 2 z cos 9 — s sgn(fc — k )^J (a 2 + a^)k 2 + 2k 2 z k{a y cosQ — a z sin#sin<^) + 



k z 



(10) 



and the spin expectation values are given by 



As can be seen from Eq. the necessary condition 



(<r x + ia y )k,s 



m 



As) 



As) 



'y}k,s 



'z/k,s 



= 1 holds for all values of 



(11) 
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FIG. 2: (Color online) Dispersion and spin expectation values on the (s = +)-branch for a magnetic field in a) x- , b) y- and 
c) z-direction, a y /a — —0.6, a z /a — —0.8, kz/ot = 0.5. The spin on the (s = — )-branch points in the opposite direction, 
i.e. {<Ji)k. s = —{&i)k,-s- The shape of the dispersion and the fc-value at which the "energy gap" becomes smallest clearly 
depends on the direction of the magnetic field. 



and k. The existence of the confinement in y-direction 
(represented by a y ) leads to a rotation of the spin out 
of the cc-y-plain into the z-direction. This indicates that 
the ratio of a y and a z is crucial for the spin direction. 

The energy dispersion Eq. (jSJ and the spin expecta- 
tion values on the (s = +)-branch are shown in Fig. [2] as 

a function of k, with k given in units of a = «/ + a 2 

and the energy in units of E a = h 2 a 2 /2m. For \k\ > a 
the spin expectation values reach their asymptotic, k- 
independent values. The spin on the (s = — )-branch 
points in the opposite direction, i.e. (<Ji}k. s = —(&i)k,-s, 
and is not shown explicitly here. In combination with 
the fact that for B = (.6,0,0), (<J y )k,s and {cr z )k,s are 
symmetric with respect to k = on both branches, this 
explains why there is no ground state magnetization in 
the y- and z-direction for B being parallel to the wire. 
However, there is a nonvanishing ground state magne- 
tization in the x-direction. The "energy gap" is given 
by 4yC(fco) [see Eq.©] and does not necessarily de- 
crease from its maximum value 4k z , if B is tilted against 
e x as stated in Ref. [9„ In units of the Zeeman energy 
Ez = 2h 2 k z /2m, the size of the "gap" Eg for arbitrary 
magnetic field B = B{sa\ 9 cos tj>, sin 9 sin </>, cos 9) is given 
by 



Eg 
Ez 



1 



a 7 sin v sin < 



(12) 



Therefore, a finite a y term is necessary for opening the 
"gap" for B||e w . To emphasize this effect, we choose 
the parameter set (a y ,a z ,kz)/a = (—0.6,-0.8,0.5) in 
Fig. [2] In many experimental systems the confining po- 
tential in the y-direction might be much weaker than 
in the z-direction. In this case \a y \ <C \a z \ but sub- 
band mixing becomes relevant. The latter strongly af- 
fects the spin-dependent transport properties as e.g. in- 
vestigated in Ref. H, and the polarization effects discussed 
here can be expected to disappear. To achieve spin po- 
larization in the present setup a strong confinement in 
the y-direction leading to a sizable a y is thus essential. 
The lower dispersion branch in Fig. [2] has a "W"-like 
shape. For B = (B,0, 0), the condition for this behavior 
is a 2 + a 2 > 2k 2 z and becomes much more complex for 



arbitrary magnetic field. We will focus on the situation 
where B = (5,0,0). 

The transmissions t sg > (conductance divided by e 2 /h) 
of an electron current at fixed Fermi energy Ep passing a 
potential step in the wire direction [see Fig. [J a)] are ob- 
tained by assuming continuity of the wave functions and 
their derivatives at the interface. Here the first index la- 
bels the branch to the left and the second index labels the 
branch to the right of the potential step. It was argued 
in Ref. [l5j that one has to consider the continuity of the 
wave function's flux and not simply its derivative, but 
in our setup both conditions lead to the same equations 
as we consider a homogeneous SOI. The total transmis- 
sion T is the sum of the four components t++, , t \-, 

and t To the right of the potential and for momenta 

|fc| > a, one can assign spins with quantum numbers |, j 
and a properly chosen quantization axis to the branches 
s = +, — because of the independence of (cr)k,s on fc. 
However, the polarization vector is given by 



t+++t- 
T 



t+-+t. 
T 



\<r)h,- 



(13) 



Since the potential step geometry was already 
discussed, 9 we will only shortly mention the influence of 
the additional term H' Rl defined in Eq. (H|), and discuss 
the interesting case of a potential barrier [see Fig.[T]b)] in 
more detail. The latter can experimentally be achieved 
by adding gates to the ID quantum wire. 

As shown in Fig. O the total polarization P = |P| of 
the current passing the potential step is large for ener- 
gies in the "gap" and increases with a. Similar to the 
transmissions t ss > , P as well as the parallel polarization 
P x depend only on Vq, kz, and a for BUe^ and not on a y 
and a z independently. The relevant energy scale of the 
polarization shown in Fig.[3]is given by Ez-, which defines 
the size of the "gap" [see Eq. (fT2)) ]. Therefore, energies 
are given in units of Ez and wave vectors in units of kz ■ 
The same holds for the transmissions and polarizations 
shown further down (see Figs. 0]and ©. The parameters 
in Fig. M are V /E z = 15, and a/k z = 2, 2.5, 3, 5. The 
energy offset is chosen such that EpjEz = corresponds 
to the middle of the "gap" . The parallel polarization P x 
gives the main contribution to the total polarization as 




FIG. 3: (Color online) Polarization of the transmission cur- 
rent at a potential step as a function of the Fermi energy for 
Vo/Ez = 15 and a/kz = 2, 2.5, 3, 5. The total polarization 
P is sizable for energies in the "gap" (indicated by the ar- 
rows). In this regime it is mostly carried by P y and P z . The 
polarization becomes negligible for energies outside the "gap" 
where P x dominates. 



the energy departs from the "gap", P x /P — > 1. How- 
ever, in this region the total polarization is negligible 
and within the "gap", the parallel component plays an 
inferior role. The ratio of the two perpendicular polar- 
izations is given by \P Z /P y \ = a> y /a z . Therefore, the 
orthogonal polarization = (0,P y ,P z ) can be rotated 
within the y-z-plane by adjusting a y and a z . 

We next study the transmission current at a potential 
barrier of height Vq and width 2x c [see Fig. [T]b)]. This 
situation might be more realistic than a simple poten- 
tial step if one thinks of further structuring by applying 
gates to the quantum wire. Fig. [4] shows the four com- 
ponents of the transmission as a function of Ep / Ez for 
a/kz = 2, 2.5, 3, Vq/Ez = 15, and kzx c = 1. Again, the 
SOI affects the transmissions t ss > only via a. Interest- 
ingly and in contrast to the potential step, the s-flipping 

transmissions are degenerate, = t |_. This can be 

understood, if one considers the possible s-flips at the 
two interfaces leading to an overall s-fiip. Labeling the 
left interface (1) and the right (2), one simply has to take 
the sum of the products of transmissions at each interface 
and obtains 



= t- 
= t. 



_(!)*_ 



-(2)+*h 



.(l)t. 

-(1)*H 



(2) 
-(2) 



(14) 



An analysis of the potential step problem shows that 
the s-conserving transmissions t++ and t are indepen- 
dent of the sign of Vo and the s-flipping transmissions 

just swap, i.e. t+- (1) = t |_(2) and t |_(1) = t-\ (2). 

This leads to exactly the same values of and t |_ 

in Eq. (fl~4|) . The exponential suppression of i++(l) and 

t |_(1) for energies within the "gap" does not affect this 

behavior. The s-conserving transmissions t ++ and t 

show an oscillatory behavior, which is well known from 



FIG. 4: (Color online) Partial transmissions at a potential 
barrier as a function of the energy for Vo/Ez = 15 and 
a/kz = 2, 2.5, 3. 
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FIG. 5: (Color online) Polarization of the transmission cur- 
rent at a potential barrier as a function of the energy for the 
same parameters as in Fig. [4] The polarization is sizable for 
energies well beyond the "gap" (indicated by the arrows) and 
shows oscillatory behavior. The ^-component P x is only rel- 
evant in regimes where the total polarization is small. 



scattering off a potential step at vanishing SOI. How- 
ever, especially for low energies, the amplitude strongly 

depends on a. The s-flipping transmissions t-\ and t |_ 

oscillate as well. The second peak of which lies in 

the "energy gap" , is suppressed compared to t , since 

right-moving (s = +)-waves are exponentially damped 
in the barrier region and therefore, as shown in Ref. 0, 

t is the dominant component at each interface in this 

energy range. 

Fig.[5]shows P and P x /P for the same parameters as in 
Fig. [U and a/kz = 2, 2.5, 3. Similarly to the potential 
step case, P = |P| and P x only depend on a and not on 
a y and a z independently. Surprisingly, the polarization 
now has a sizable value in an energy interval much bigger 
than the "gap" , which just goes from — Ez to Ez (see the 
arrows in Fig. [S]) . This behavior must be contrasted to 
the polarization in the case of a potential step as shown in 
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Fig.[3]and first introduced in Ref.[§. It can be traced back 

to the energy dependence of i H and t |_ shown in Fig. 2] 

Both have finite weight well beyond the "energy gap". 
This might be due to interference effects of transmitted 
and reflected waves in the barrier region. 



III. LATTICE MODEL 

In a next step, we are aiming at constructing a tight- 
binding lattice model which in appropriate parameter 
regimes shows similar physics as our continuum model. 
This will put us in a position to study the effect of 
electron-electron interaction neglected so far using the 
functional renormalization group method^ In ID wires 
the two-particle interaction is known to strongly alter 
the low-energy physics of many-body systems leading to 
so called Luttinger liquid behavior^ It can be expected 
that the interplay of the SOI effects discussed above and 
correlation effects leads to interesting physics. The SOI 
can be modeled by spin-flip hopping terms with ampli- 
tude a y and a z in a usual tight-binding models 

We start with a representation of the Hamiltonian in 
terms of Wannier states \j, a) with j £ Z labeling the 
lattice site and a =f , J. labeling the spin. The spin quan- 
tization is chosen along the z-direction. With ct a being 
the creation operator of an electron at site j with spin er, 
the lattice model Hamiltonian for an arbitrary magnetic 
field B = B (sin 9 cos ip, sin 9 sin ip, cos 9) can be written 
as 



H — H o + Hpot + H R + H z , 



(15) 



with the free part 

(16) 

containing the on-site energy and the conventional (spin- 
conserving) hopping, external potential (due to e.g. nano- 
device structuring) 



-Hpot = Yl ^>4,< 



•Cj.cr , 



(17) 



the spin-flip (Rashba) hopping terms 

H R = -a z Y (c]+i,<t ( i<7 v)a,*> C J>' +H.c.) (18) 



and the Zeeman term 



Hz = 2k\ Y C L ( a *)*,*' sin6 cost P 



(19) 



+ ( a y)a,<r> sin 6> sin 93 + {<J z ) aa , cosfrj c,>/ 
We show the analogy to the continuum case suppressing 
H pot and take as an ansatz for the corresponding eigen- 
states 



\k,s)=J2<(k)e ikj \j,*) 



(20) 



This leads to the eigenenergies 



E {s) (k) = e- It cos k + 2s sgn(fc - k )^D(k) , (21) 



with 



k 



and 



-k%{a y t 



a z sm 



>smcp)/{a 2 y + a 2 z )] 



(22) 



D(k) 



(a 2 , + a 2 z ) sin 2 k + k% (23) 
+2k z sin k (a y cos 9 — a z sin 9 sin (p) . 



Eq. (|2~lj) has almost the same form as the continuum ver- 
sion Eq. ©. In fact, choosing the on-site energy e = 2t, 
which corresponds just to an overall energy shift, and 
substituting cosfc by 1 — fc 2 /2 and sinfc by k, which is 
valid for sufficiently small |fc|, we get exactly the same 
form. Note however that, in contrast to the continuum 
case, a y ,a z and k 2 z now have the unit of energy, but 

since ct , defined in Eq. is dimensionless, all formu- 
las remain valid. We choose for the eigenstates Eq. (|2"0|) 
af(fc) = 1 and obtain a|(fc) = cj^ 



with 



» 



-ia z sin k 



k\e 



■up , 



sin/c + k% cos 9 — s sgn(fc — k )J (a y + a 2 z ) sin k + 2k 2 z sin k(a y cos 9 — a z sinflsin^) + k z 



(24) 
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FIG. 6: (Color online) Lattice dispersion and spin expectation values on the (s = +)-branch for a magnetic field in a) x- , b) y- 
and c) 2-direction for t/a = 1, a y /a — —0.6, a z /a — —0.8, kz/a — 0.5. The spin on the (s = — )-branch points in the opposite 
direction, i.e. (<7i)fc ]S = — {(?i}k,-s- For |fe — feo| < t/2 one obtains exactly the same behavior as in the continuum case. 



and the spin expectation values have exactly the contin- 
uum form 



(<7 X + i<Jy)k,s = 2- 



-1 + 



1 + 



\Vz)k.s = 



» 



1 



» 



The energy dispersions and the spin expectation values 
for magnetic fields in x-, y-, and z-direction are shown 
in Fig. [6] Besides the cosine-like structure, which be- 
comes especially relevant near the band edges, the dis- 
persion and spin expectation values have the same shape 
as in the continuum model. A direct comparison of Fig. [5] 
and Fig. [2] shows that our lattice model reproduces the 
low energy physics, i.e. for \k — ko\ < 7r/2, observed in 
the continuum. As above we only show the spin ex- 
pectation values on the (s = +)-branch. The spin on 
the (s = — )-branch points in the opposite direction, 
i.e. (<Ji)k,s — ~ (fi)fc.-s- The direct relation between the 
dispersion and the spin expectation values for energies 
of the order of the "gap" is the essential feature leading 
to the remarkable scattering properties of the continuum 
model (and eventually a spin polarized conductance) at 
steps and barriers. One can thus expect similar trans- 
port characteristics to be realized in the lattice model. A 
detailed discussion of this and in particular the effect of 
the electron-electron interaction on transport will be the 
topic of an upcoming publication^ 



as arbitrary magnetic field, and have shown that an ad- 
ditional SOI term resulting from the lateral confinement 
of a 2DEG to a ID wire geometry leads to a rotation 
of the spin out of the 2D plane. For the case of a mag- 
netic field parallel to the quantum wire, the transmission 
and polarization of a linear response current at a poten- 
tial step as well as at a potential barrier were studied. 
For the latter, we observed an extended energy range, 
where significant spin polarization can be achieved. We 
showed that this spin polarization can be rotated out of 
the plane of the 2DEG arbitrarily by adjusting the SOI 
constants a y and a z . The potential barrier describes a 
setup which can experimentally be achieved by adding 
further gates to the wire geometry. We then constructed 
a lattice model which shows the same low energy physics 
as the continuum model. This lattice model now enables 
us to investigate the interplay of SOI and Coulomb inter- 
action in quantum wires with potential steps and barriers 
using the functional renormalization group method i 13 ' 14 
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